Measurement of CP-violation asymmetries in D° — >■ 7r+ n 

., 17 J. Antos, 12 

15 T. Arisawa, 54 A. Artikov, 13 J. Asaadi, 49 W. Ashmanskas, 15 B. Auerbach 

. 23 



G 



mv", J. Bellinger, 

yamm," A. Beretvas, i,J A. Bhatti,' 1 " U. Biselky 1 ,™ 1. tSizjaV" K.K. tSland^ B. Blumenfeld, 23 A. Bocci, 14 
Bodek, 45 D. Bortoletto, 44 J. Boudreau, 43 A. Boveia, 11 L. Brigliadori 66 , 6 C. Bromberg, 33 E. Brucken, 21 
igov, 13 H.S. Budd, 45 K. Burkett, 15 G. Busetto^, 40 P. Bussey, 19 A. Buzatu, 31 A. Calamba, 10 C. Calancha, 29 
S. Camarda, 4 M. Campanelli, 28 M. Campbell, 32 F. Canelli, 11 ' 15 B. Carls, 22 D. Carlsmith, 56 R. Carosi, 42 
. Carrillo™, 16 S. Carron, 15 B. Casal fc , 9 M. Casarsa, 50 A. Castro 66 , 6 P. Catastini, 20 D. Cauz, 50 V. Cavaliere, 22 

M. Cavalli-Sforza, 4 A. Cerri^, 26 L. Cerrito 5 , 28 Y.C. Chen, 1 M. Chertok, 7 G. Chiarelli, 42 G. Chlachidze, 15 
'. Chlebana, 15 K. Cho, 25 D. Chokheli, 13 W.H. Chung, 56 Y.S. Chung, 45 M.A. Ciocci' 1 ' 1 , 42 A. Clark, 18 C. Clarke, 55 
Compostella^, 40 M.E. Convery, 15 J. Conway, 7 M.Corbo, 15 M. Cordelli, 17 C.A. Cox, 7 D.J. Cox, 7 F. Crescioli 99 , 42 
F. Cuevas*, 9 R. Culbertson, 15 D. Dagenhart, 15 N. d'Ascenzo 1 ", 15 M. Datta, 15 P. de Barbara, 45 M. Dell'Orso 99 , 42 
L. Dcmortier, 46 M. Deninno, 6 F. Devoto, 21 M. d'Errico //40 A. Di Canto 99 , 42 B. Di Ruzza, 15 J.R. Dittmann, 5 
M. D'Onofrio, 27 S. Donati 99 , 42 P. Dong, 15 M. Dorigo, 50 T. Dorigo, 40 K. Ebina, 54 A. Elagin, 49 A. Eppig, 32 
R. Erbacher, 7 S. Errede, 22 N. Ershaidat dd , 15 R. Eusebi, 49 S. Farrington, 39 M. Feindt, 24 J.P. Fernandez, 29 
R. Field, 16 G. Flanagan", 15 R. Forrest, 7 M.J. Frank, 5 M. Franklin, 20 J.C. Freeman, 15 Y. Funakoshi, 54 I. Furic, 16 

M. Gallinaro, 46 J.E. Garcia, 18 A.F. Garfmkel, 44 P. Garosi' 1 ' 1 , 42 H. Gerberich, 22 E. Gerchtein, 15 S. Giagu, 47 
V. Giakoumopoulou, 3 P. Giannetti, 42 K. Gibson, 43 CM. Ginsburg, 15 N. Giokaris, 3 P. Giromini, 17 G. Giurgiu, 23 
V. Glagolev, 13 D. Glenzinski, 15 M. Gold, 35 D. Goldin, 49 N. Goldschmidt, 16 A. Golossanov, 15 G. Gomez, 9 
G. Gomez-Ceballos, 30 M. Goncharov, 30 O. Gonzalez, 29 I. Gorelov, 35 A.T. Goshaw, 14 K. Goulianos, 46 S. Grinstein, 4 

C. Grosso-Pilcher, 11 R.C. Group 53 , 15 J. Guimaraes da Costa, 20 S.R. Hahn, 15 E. Halkiadakis, 48 A. Hamaguchi, 38 
JY. Han, 45 F. Happacher, 17 K. Hara, 51 D. Hare, 48 M. Hare, 52 R.F. Harr, 55 K. Hatakeyama, 5 C. Hays, 39 M. Heck, 24 
J. Heinrich, 41 M. Herndon, 56 S. Hewamanage, 5 A. Hocker, 15 W. Hopkins 9 , 15 D. Horn, 24 S. Hou, 1 R.E. Hughes, 36 
'.. Hurwitz, 11 U. Husemann, 57 N. Hussain, 31 M. Hussein, 33 J. Huston, 33 G. Introzzi, 42 M. Iori", 47 A. Ivanov p , 7 

E. James, 15 D. Jang, 10 B. Jayatilaka, 14 E.J. Jeon, 25 S. Jindariani, 15 M. Jones, 44 K.K. Joo, 25 SY. Jun, 10 
I Junk, 15 T. Kamon 25 , 49 P.E. Karchin, 55 A. Kasmi, 5 Y. Kato°, 38 W. Ketchum, 11 J. Keung, 41 V. Khotilovich, 49 
B. Kilminster, 15 D.H. Kim, 25 H.S. Kim, 25 J.E. Kim, 25 M.J. Kim, 17 S.B. Kim, 25 S.H. Kim, 51 Y.K. Kim, 11 
x.J. Kim, 25 N. Kimura, 54 M. Kirby, 15 S. Klimenko, 16 K. Knoepfel, 15 K. Kondo*, 54 D.J. Kong, 25 J. Konigsberg, 16 

A.V. Kotwal, 14 M. Kreps mm , 24 J. Kroll, 41 D. Krop, 11 M. Kruse, 14 V. Krutelyov 6 , 49 T. Kuhr, 24 M. Kurata, 51 
S. Kwang, 11 A.T. Laasanen, 44 S. Lami, 42 S. Lammel, 15 M. Lancaster, 28 R.L. Lander, 7 K. Lannon 9 , 36 A. Lath, 48 
G. Latino' 1 ' 1 , 42 T. LeCompte, 2 E. Lee, 49 H.S. Lee*, 11 J.S. Lee, 25 S.W. Lee hb , 49 S. Leo 99 , 42 S. Leone, 42 J.D. Lewis, 15 
A. Limosani*, 14 C.-J. Lin, 26 M. Lindgren, 15 E. Lipeles, 41 A. Lister, 18 D.O. Litvintsev, 15 C. Liu, 43 H. Liu, 53 Q. Liu, 44 
T. Liu, 15 S. Lockwitz, 57 A. Loginov, 57 D. Lucchesi // , 40 J. Lueck, 24 P. Lujan, 26 P. Lukens, 15 G. Lungu, 46 J. Lys, 26 
R. Lysak 6 , 12 R. Madrak, 15 K. Maeshima, 15 P. Maestro' 1 ' 1 , 42 S. Malik, 46 G. Manca a , 27 A. Manousakis-Katsikakis, 3 
F. Margaroli, 47 C. Marino, 24 M. Martinez, 4 P. Mastrandrea, 47 K. Matera, 22 M.E. Mattson, 55 A. Mazzacane, 15 

P. Mazzanti, 6 K.S. McFarland, 45 P. Mclntyrc, 49 R. McNulty-?, 27 A. Mehta, 27 P. Mehtala, 21 C. Mesropian, 46 
T. Miao, 15 D. Mietlicki, 32 A. Mitra, 1 H. Miyake, 51 S. Moed, 15 N. Moggi, 6 M.N. Mondragon™, 15 C.S. Moon, 25 
R. Moore, 15 M.J. Morello", 42 J. Morlock, 24 P. Movilla Fernandez, 15 A. Mukherjee, 15 Th. Muller, 24 P. Murat, 15 

M. Mussini ee , 6 J. Nachtman™, 15 Y. Nagai, 51 J. Naganoma, 54 I. Nakano, 37 A. Napier, 52 J. Nett, 49 C. Neu, 53 
M.S. Neubauer, 22 J. Nielsen^, 26 L. Nodulman, 2 SY. Noh, 25 O. Norniella, 22 L. Oakes, 39 S.H. Oh, 14 Y.D. Oh, 25 
I. Oksuzian, 53 T. Okusawa, 38 R. Orava, 21 L. Ortolan, 4 S. Pagan Griso // , 40 C. Pagliarone, 50 E. Palencia / , 9 
V. Papadimitriou, 15 A. A. Paramonov, 2 J. Patrick, 15 G. Pauletta fcfc , 50 M. Paulini, 10 C. Paus, 30 D.E. Pellett, 7 
A. Penzo, 50 T.J. Phillips, 14 G. Piacentino, 42 E. Pianori, 41 J. Pilot, 36 K. Pitts, 22 C. Plager, 8 L. Pondrom, 56 
S. Poprocki 9 , 15 K. Potamianos, 44 F. Prokoshin cc , 13 A. Pranko, 26 F. Ptohos' 1 , 17 G. Punzi 99 , 42 A. Rahaman, 43 

V. Ramakrishnan, 56 N. Ranjan, 44 I. Redondo, 29 P. Renton, 39 M. Rescigno, 47 T. Riddick, 28 F. Rimondi ee , 6 
L. Ristori 42 , 15 A. Robson, 19 T. Rodrigo, 9 T. Rodriguez, 41 E. Rogers, 22 S. Rolli*, 52 R. Roser, 15 F. Ruffini' 1 ' 1 , 42 
A. Ruiz, 9 J. Russ, 10 V. Rusu, 15 A. Safonov, 49 W.K. Sakumoto, 45 Y. Sakurai, 54 L. Santi fefe , 50 K. Sato, 51 
V. Saveliev'V 5 A. Savoy-Navarro aa , 15 P. Schlabach, 15 A. Schmidt, 24 E.E. Schmidt, 15 T. Schwarz, 15 L. Scodellaro, 9 
A. Scribano' 1 ' 1 , 42 F. Scuri, 42 S. Seidel, 35 Y. Seiya, 38 A. Semenov, 13 F. Sforza' 1 ' 1 , 42 S.Z. Shalhout, 7 T. Shears, 27 
P.F. Shepard, 43 M. Shimojima", 51 M. Shochet, 11 I. Shreyber-Tecker, 34 A. Simonenko, 13 P. Sinervo, 31 K. Sliwa, 52 



JY 
J. 

M 

T.R 
Y.J 



2 



J.R. Smith, 7 F.D. Snider, 15 A. Soha, 15 V. Sorin, 4 H. Song, 43 P. Squillacioti' 1 ' 1 , 42 M. Stancari, 15 R. St. Denis, 19 
B. Stelzer, 31 O. Stclzer-Chilton, 31 D. Stentz*, 15 J. Strologas, 35 G.L. Strycker, 32 Y. Sudo, 51 A. Sukhanov, 15 
I. Suslov, 13 K. Takemasa, 51 Y. Takeuchi, 51 J. Tang, 11 M. Tecchio, 32 P.K. Teng, 1 J. Thom ff , 15 J. Thome, 10 
G.A. Thompson, 22 E. Thomson, 41 D. Toback, 49 S. Tokar, 12 K. Tollefson, 33 T. Tomura, 51 D. Tonelli, 15 S. Torre, 17 
D. Torretta, 15 P. Totaro, 40 M. Trovato", 42 F. Ukegawa, 51 S. Uozumi, 25 A. Varganov, 32 F. Vazquez" 1 , 16 G. Velev, 15 
C. Vellidis, 15 M. Vidal, 44 I. Vila, 9 R. Vilar, 9 J. Vizan, 9 M. Vogel, 35 G. Volpi, 17 P. Wagner, 41 R.L. Wagner, 15 
T. Wakisaka, 38 R. Wallny, 8 S.M. Wang, 1 A. Warburton, 31 D. Waters, 28 W.C. Wester III, 15 D. Whiteson 6 , 41 
A.B. Wicklund, 2 E. Wickhmd, 15 S. Wilbur, 11 F. Wick, 24 H.H. Williams, 41 J.S. Wilson, 36 P. Wilson, 15 
B.L. Winer, 36 P. WittichV 5 S. Wolbers, 15 H. Wolfe, 36 T. Wright, 32 X. Wu, 18 Z. Wu, 5 K. Yamamoto, 38 
D. Yamato, 38 T. Yang, 15 U.K. Yang r , n Y.C. Yang, 25 W.-M. Yao, 26 G.P. Ych, 15 K. Yi n , 15 J. Yoh, 15 K. Yorita, 54 
T. Yoshida', 38 G.B. Yu, 14 I. Yu, 25 S.S. Yu, 15 J.C. Yun, 15 A. Zanetti, 50 Y. Zcng, 14 C. Zhou, 14 and S. Zucchelli ee6 

(CDF Collaboration^) 

1 Institute of Physics, Accidentia Sinica, Taipei, Taiwan 11529, Republic of China 
2 Argonne National Laboratory, Argonne, Illinois 60439, USA 
s University of Athens, 157 71 Athens, Greece 
4 Institut de Fisica d'Altes Energies, ICREA, Universitat Autonoma de Barcelona, E-08193, Bellaterra (Barcelona), Spain 

5 Baylor University, Waco, Texas 76798, USA 
6 Istituto Nazionale di Fisica Nucleare Bologna, ee University of Bologna, 1-40127 Bologna, Italy 
7 University of California, Davis, Davis, California 95616, USA 
8 University of California, Los Angeles, Los Angeles, California 90024, USA 
9 Institute de Fisica de Cantabria, CSIC-University of Cantabria, 39005 Santander, Spam 
10 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
11 Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637, USA 
12 Cornenius University, 842 48 Bratislava, Slovakia; Institute of Experimental Physics, 040 01 Kosice, Slovakia 
13 Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 

14 Duke University, Durham, North Carolina 27708, USA 
Fermi National Accelerator Laboratory, Batavia, Illinois 60510, USA 
16 University of Florida, Gainesville, Florida 32611, USA 
17 Laboratori Nazionali di Frascati, Istituto Nazionale di Fisica Nucleare, 1-00044 Frascati, Italy 
18 University of Geneva, CH-1211 Geneva 4, Switzerland 
19 Glasgow University, Glasgow G12 8QQ, United Kingdom 
20 Harvard University, Cambridge, Massachusetts 02138, USA 
Division of High Energy Physics, Department of Physics, 
University of Helsinki and Helsinki Institute of Physics, FIN-00014, Helsinki, Finland 
22 University of Illinois, Urbana, Illinois 61801, USA 
28 The Johns Hopkins University, Baltimore, Maryland 21218, USA 
2J * Institut fur Experimentelle Kernphysik, Karlsruhe Institute of Technology, D-76131 Karlsruhe, Germany 
25 Center for High Energy Physics: Kyungpook National University, 
Daegu 702-701, Korea; Seoul National University, Seoul 151-742, 
Korea; Sungkyunkwan University, Suwon 440-746, 
Korea; Korea Institute of Science and Technology Information, 
Daejeon 305-806, Korea; Chonnam National University, Gwangju 500-757, 
Korea; Chonbuk National University, Jeonju 561-756, Korea 
26 Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 
27 University of Liverpool, Liverpool L69 7ZE, United Kingdom 
28 University College London, London WC1E 6BT, United Kingdom 
29 Centro de Investigaciones Energeticas Medioambientales y Tecnologicas, E-28040 Madrid, Spain 
80 Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 
31 Institute of Particle Physics: McGill University, Montreal, Quebec, 
Canada H3A 2T8; Simon Fraser University, Burnaby, British Columbia, 
Canada V5A 1S6; University of Toronto, Toronto, Ontario, 
Canada M5S 1A7; and TRIUMF, Vancouver, British Columbia, Canada V6T 2A3 
32 University of Michigan, Ann Arbor, Michigan 48109, USA 
33 Michigan State University, East Lansing, Michigan 48824, USA 
34 Institution for Theoretical and Experimental Physics, ITEP, Moscow 117259, Russia 
35 University of New Mexico, Albuquerque, New Mexico 87131, USA 
36 The Ohio State University, Columbus, Ohio 43210, USA 
37 Okayama University, Okayama 700-8530, Japan 
38 Osaka City University, Osaka 588, Japan 
39 University of Oxford, Oxford 0X1 3RH, United Kingdom 
40 Istituto Nazionale di Fisica Nucleare, Sezione di Padova-Trento, ^University of Padova, 1-35131 Padova, Italy 



3 



41 University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 

42 Istituto Nazionale di Fisica Nucleare Pisa, 99 University of Pisa, 
University of Siena and "Scuola Normale Superiore, 1-56127 Pisa, Italy 
4 3 University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 
44 Purdue University, West Lafayette, Indiana 47907, USA 
45 University of Rochester, Rochester, New York 14-627, USA 
46 The Rockefeller University, New York, New York 10065, USA 
4 'istituto Nazionale di Fisica Nucleare, Sezione di Roma 1, 

33 Sapienza Universita di Roma, 1-00185 Roma, Italy 
48 Rutgers University, Piscataway, New Jersey 08855, USA 
49 Texas A&M University, College Station, Texas 77843, USA 

50 Istituto Nazionale di Fisica Nucleare Trieste/ U dine, 
1-34100 Trieste, kk University of Udme, 1-33100 Udme, Italy 
51 University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
52 Tufts University, Medford, Massachusetts 02155, USA 
53 University of Virginia, Charlottesville, Virginia 22906, USA 
54 Waseda University, Tokyo 169, Japan 
Wayne State University, Detroit, Michigan 48201, USA 
56 University of Wisconsin, Madison, Wisconsin 53706, USA 
Yale University, New Haven, Connecticut 06520, USA 

We report a measurement of time-integrated CP-violation asymmetries in the resonant substruc- 
ture of the three-body decay D° — > Kg-K + ir~ using CDF II data corresponding to 6.0 fb -1 of 
integrated luminosity from Tevatron pp collisions at y/s = 1.96 TeV. The charm mesons used in 
this analysis come from D* + (2010) — > D°tt + and D* - (2010) — ¥ D°tv~ , where the production flavor 
of the charm meson is determined by the charge of the accompanying pion. We apply a Dalitz- 
amplitude analysis for the description of the dynamic decay structure and use two complementary 
approaches, namely a full Dalitz-plot fit employing the isobar model for the contributing resonances 
and a model-independent bin-by-bin comparison of the D° and D° Dalitz plots. We find no CP- 
violation effects and measure an asymmetry of Acp = (—0.05 ± 0.57(stat) ± 0.54(syst))% for the 
overall integrated CP-violation asymmetry, consistent with the standard model prediction. 

PACS numbers: 13. 25. Ft, 11.30.Er 



I. INTRODUCTION 



'Deceased 

tWith visitors from "Istituto Nazionale di Fisica Nucleare, Sczionc 
di Cagliari, 09042 Monserrato (Cagliari), Italy, ^University of CA 
Irvine, Irvine, CA 92697, USA, c University of CA Santa Barbara, 
Santa Barbara, CA 93106, USA, d University of CA Santa Cruz, 
Santa Cruz, CA 95064, USA, institute of Physics, Academy of 
Sciences of the Czech Republic, Czech Republic, ^CERN, CH-1211 
Geneva, Switzerland, ^Cornell University, Ithaca, NY 14853, USA, 
^University of Cyprus, Nicosia CY-1678, Cyprus, 'Office of Sci- 
ence, U.S. Department of Energy, Washington, DC 20585, USA, 
J University College Dublin, Dublin 4, Ireland, fe ETH, 8092 Zurich, 
Switzerland, 'University of Fukui, Fukui City, Fukui Prefecture, 
Japan 910-0017, m Universidad Iberoamericana, Mexico D.F., Mex- 
ico, "University of Iowa, Iowa City, IA 52242, USA, °Kinki Uni- 
versity, Higashi-Osaka City, Japan 577-8502, ^Kansas State Uni- 
versity, Manhattan, KS 66506, USA, 9 Ewha Womans University, 
Seoul, 120-750, Korea, ""University of Manchester, Manchester M13 
9PL, United Kingdom, s Queen Mary, University of London, Lon- 
don, El 4NS, United Kingdom, *University of Melbourne, Victoria 
3010, Australia, "Muons, Inc., Batavia, IL 60510, USA, "Nagasaki 
Institute of Applied Science, Nagasaki, Japan, w National Research 
Nuclear University, Moscow, Russia, x Northwestern University, 
Evanston, IL 60208, USA, ^University of Notre Dame, Notre Dame, 
IN 46556, USA, z Universidad de Oviedo, E-33007 Ovicdo, Spain, 
aa CNRS-IN2P3, Paris, F-75205 France, 6i, Texas Tech University, 
Lubbock, TX 79609, USA, cc Universidad Tecnica Federico Santa 
Maria, llOv Valparaiso, Chile, dti Yarmouk University, Irbid 211- 
63, Jordan, mm University of Warwick, Coventry CV4 7AL, United 



The phenomenon of CP violation is well established for 
weakly decaying hadrons consisting of down-type quarks. 
It has been observed with mixing (indirect CP violation) 
and without mixing (direct CP violation) in decays of the 
K° and B Q mesons, and confirms the theory of Kobayashi 
and Maskawa pj that describes CP violation in the stan- 
dard model. For charm mesons, CP-violation effects are 
expected to be small. The lack of experimental observa- 
tion of CP violation in charm meson decays is consistent 
with the expectation. Only recently the LHCb collabo- 
ration has reported the first evidence for CP violation 
in the time-integrated rate of D° decays to two hadrons 
atO(10" 2 ) 0, which was consequently confirmed by the 
CDF collaboration [jjj]. Whether this is consistent with 
the standard model expectation or a hint for new physics 
is not yet clear. Therefore, it is important to complement 
the LHCb result with measurements of CP asymmetries 
in other D° decay modes. 

In this article, we describe a search for CP viola- 
tion in time-integrated decay rates of D° mesons to the 
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Kg 7r + 7r~ final state. The standard model expectations 
for the size of time-integrated CP- violation asymmetries 
in Kgir + ir~ decays are 0(1O~ 6 ), where the dominant 
contribution arises from CP violation in K°-K° mix- 
ing 0, @. The observation of any significantly larger 
asymmetries would be a strong hint for physics beyond 
the standard model. 

The CLEO collaboration performed a dedicated search 
for time-integrated CP violation in a fit to the Dalitz plot 
of the D°/D° three-body decay to K% tt+ tt" ||. Bcnc 
and BABAR allowed for CP violation in their measure- 
ments of the D°-D° mixing parameters @J1| ■ Up to now, 
no CP-violation effects have been found [0). 

In this analysis, we exploit a large sample of D* (2010) ± 
decays, reconstructed in a dataset corresponding to 
6.0 fb -1 of integrated luminosity produced in pp col- 
lisions at yfs = 1.96 TeV and collected by the CDF 
II detector, to measure time-integrated CP- violation 
asymmetries in the resonant substructure of the decay 
D°/D° — > Kg7r + 7r~. The neutral D meson produc- 
tion flavor is determined by the charge of the pion in 
the D*+(2010) D°tt+ and P>*-(2010) -» D ^- de- 
cays (D* tagging). Throughout the rest of the paper 
the use of a specific particle state implies the use of the 
charge-conjugate state as well, unless explicitly noted. 
For brevity, D*(2010)+ is abbreviated as D*+ . 

In Sections |TT] and IIIII we briefly describe the CDF II 
detector and the trigger components important for this 
analysis. In Sections IIVI and [V] we describe the offline 
candidate reconstruction and selection, respectively. In 
Sec. I VII we explain the Dalitz-plot fits of the resonant 
substructure of the decay, followed by a discussion of sys- 
tematic uncertainties in Sec. I VIII and a presentation of 
the results in Sec. I Villi In Sec. |IX] we describe a model- 
independent search for CP asymmetries using a bin-by- 
bin comparison of Dalitz plots [l(| , followed by the con- 
clusion in Sec. [X] 



II. CDF II DETECTOR 

The analysis is performed on a dataset collected with 
the CDF II detector [TO] between February 2002 and 
February 2010, corresponding to an integrated luminosity 
of 6.0fb _1 of collisions. Among the components and 
capabilities of the detector, the charged particle tracking 
is the one most relevant to this analysis. The tracking 
system lies within a uniform, axial magnetic field of 1.4 T. 
The inner tracking volume, up to a radius of 28 cm, is 
composed of six or seven layers, depending on polar an- 
gle, of double-sided silicon microstrip detectors [l2l. [l3|. 
An additional single-sided silicon layer is mounted di- 
rectly on the beam pipe at a radius of 1.5 cm [l4j]> al- 
lowing excellent resolution of the transverse impact pa- 
rameter do, defined as the distance of closest approach 
of a charged particle trajectory (track) to the interac- 
tion point in the plane transverse to the beam line. The 
silicon detector allows the identification of displaced, sec- 



ondary vertices from bottom-hadron and charm-hadron 
decays with a resolution of approximately 30 pm in the 
transverse and 70 pm in the longitudinal direction. The 
outer tracking volume from a radius of 40 to 137 cm is 
occupied by an open-cell argon-ethane gas drift cham- 
ber (COT) |l5[. An important aspect for this analysis 
is that the layout of the COT is intrinsically charge- 
asymmetric because of an about 35° tilt angle between 
the cell orientation and the radial direction. The total 
tracking system provides a transverse momentum reso- 
lution of a{p T )/p T « 0.07% (GeV/cT 1 for tracks with 
Pt > 2 GeV/c. A more detailed description of the track- 
ing system can be found in Ref. (l6j . 



III. ONLINE EVENT SELECTION 

A three-level event-selection system (trigger) is used. 
At level 1, a hardware track-processor [17| identifies 
charged particles using information from the COT and 
measures their transverse momenta and azimuthal an- 
gles around the beam direction. The basic requirement 
at level 1 is the presence of two charged particles, each 
with px > 2 GcV/c. At level 2. the silicon vertex trig- 
ger [181 ] adds silicon-hit information to the tracks found 
by the hardware track-processor, thus allowing the pre- 
cise measurement of impact parameters of tracks. The 
two level- 1 tracks are required to have impact parame- 
ters between 0.1 and 1 mm, an opening angle in the trans- 
verse plane between 2° and 90°, and to be consistent with 
coming from a common vertex displaced from the inter- 
action point by at least 200 pm in the plane transverse 
to the beam line. This is complemented by a selection 
without the vertex displacement requirement, collecting 
events with low invariant-mass track pairs having open- 
ing angles less than 6° . The level-3 trigger is implemented 
in software and provides the final online selection by con- 
firming the first two trigger-level decisions using a more 
precise reconstruction. This trigger is designed to col- 
lect hadronic decays of long-lived particles such as b and 
c hadrons. Three different configurations of this trigger 
are employed, requiring a minimum on the scalar sum of 
the transverse momenta of the two trigger tracks of 4.0, 
5.5, and 6.5 GeV/c. The active threshold depends on the 
instantaneous luminosity conditions, with higher thresh- 
olds used at higher instantaneous luminosity to reduce 
the higher trigger accept-rate. 



IV. OFFLINE EVENT RECONSTRUCTION 

The offline reconstruction of candidates starts by fit- 
ting tracks taking into account multiple scattering and 
ionization energy-loss calculated for the pion mass hy- 
pothesis. Since all final state particles in the studied de- 
cay chain are pions, we assign the pion mass to each track 
in the following steps. Two oppositely-charged tracks 
are combined to form a Kg candidate. To construct 
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D° candidates, each K s candidate is combined with all 
other possible oppositely charged track pairs found in 
the event. Finally, the D* + candidates are obtained by 
combining each D° candidate with one of the remaining 
tracks in the event. The tracks forming the Kg, D°, and 
D* + candidates are subjected to separate kinematic fits 
that constrain them to originate from a common decay 
point in each case, resulting in three reconstructed decay 
points for the considered decay chain. In each step of the 
reconstruction, standard quality requirements on tracks 
and vertices are used to ensure well-measured momenta 
and decay positions [l9| . 



V. CANDIDATE SELECTION 

For the selection of the candidates, we first impose 
some quality requirements to suppress the most obvi- 
ous backgrounds, such as random combinations of low- 
Pt particles, D* + pions strongly displaced from the in- 
teraction point, or short-lived D° candidates without 
proper secondary-vertex separation. We require the 
transverse momentum of each pion to be greater than 
400MeV/c, the transverse momentum of the D* + to ex- 
ceed 5.0GeV/c, the impact parameter of the pion from 
the D* + decay divided by its uncertainty do/o'do ( 7r .D*+ ) 
not to exceed 15, and the transverse decay length of the 
D° candidate projected into the transverse momentum 
{L xy ) to exceed its resolution (ox ). For the surviv- 
ing candidates we use an artificial neural network to dis- 
tinguish signal from background. The neural network 
is constructed with the NeuroBayes package [2(| and 
trained, using data only, by means of the s Plot tech- 
nique [21(. This technique assigns a weight to each can- 
didate, proportional to the probability that the candidate 
is signal. The candidate weight is based on a discrimi- 
nating variable, which is required to be independent of 
the ones used in the neural network training. In our 
case, the discriminating variable is the mass difference 
AM = M{K s tt+7t~ 7T+) - M(K s tt + tt-) of the D*+ can- 
didate. In the training, each candidate enters with a 
weight calculated from the signal probability that is de- 
rived from its mass. Based on these weights, the neu- 
ral network learns the features of signal and background 
events. Since we use only data for the neural network 
training, we randomly split each sample into two parts 
evenly distributed in data taking time and train two net- 
works. Each network is then applied to the complemen- 
tary subsample in order for the selection to be trained on 
a sample independent from the one to which it is applied. 
This approach avoids a bias of the selection originating 
from statistical fluctuations possibly learned by the net- 
work. The method of NeuroBayes /Plot trainings was 
first applied in our previous work on charm baryons p2j . 

The network uses five input variables. Ordered by 
decreasing importance, these are L xy /a L:cy (D ), the \ 2 
quality of the D* + vertex fit, do/ ' (Jd a {^D*+)^ Pt(tt_d»+), 
and the reconstructed mass of the K s candidate. The 



D* + network training is based on the mass difference 
distribution in the range 140 < AM < 156MeV/c 2 . A 
fit of a non-relativistic Breit-Wigner convoluted with a 
Gaussian for the signal and a third-order polynomial for 
the background function defines the probability density 
functions used to calculate the s Vlot weights. The final 
neural network output requirement is chosen to maxi- 
mize SyV S + B where S (B) is the estimated number of 
signal (background) events in the signal region estimated 
from a fit to the M(K S ir + ir~) distribution. The selected 
requirement corresponds to an a posteriori signal prob- 
ability greater than 25%. In 11% of the selected events, 
multiple D* + candidates are reconstructed for a single 
D° candidate. We choose the D* + candidate that gives 
the most D*-like neural network output and remove all 
others to avoid identical candidates populating the Dalitz 
space. 

The M(K s Tr + n~) and AM distributions of the se- 
lected candidates are shown in Fig. [TJ together with 
the corresponding fits to determine the signal and back- 
ground yields. In the K s tt + tt~ mass distribution, the D° 
signal is described by the sum of two Gaussian functions 
with common mean, and the background is modeled by 
a linear function. The D* + signal in the AM distribu- 
tion is described by a non-relativistic Breit-Wigner con- 
voluted with a resolution function. The latter is deter- 
mined from simulated events and consists of the weighted 
sum of three Gaussian functions to model the more com- 
plicated shape of the resolution close to the kinematic 
threshold in the AM distribution. The background is 
modeled by a third-order polynomial. 

For the Dalitz-plot studies, the analysis is restricted 
to candidates populating two mass ranges, 1.84 < 
M(K s tt+tt-) < 1.89GeV/c 2 and 143.4 < AM < 
147.4 MeV/c 2 , indicated by the dashed vertical lines in 
Fig. [TJ The selected data sample contains approximately 
3.5 x 10 5 signal events and consists of about 90% correctly 
.D*-tagged D° signal, 1% mistagged D° signal, and 9% 
background candidates. A mistag results from the com- 
bination of a D° with a random, wrongly-charged pion 
not originating from a D*^ decay. The mistag fraction 
is estimated by a dedicated parameter in the Dalitz-plot 
fit described in Sec. IVI Al 

The resonant substructure of a three-body decay can 
be described by means of a Dalitz-amplitude analy- 
sis (23[. The Dalitz plot of the considered decay D° — > 
K s 7T+ 7r~, composed of all selected candidates, is shown 
in Fig. [5J It represents the decay dynamics as a func- 
tion of the squared invariant masses of the two-body 
combinations A"^7r ± (RS) and 7r + 7r~, where the notation 
A'g7r ± (RS) expresses that the Cabibbo-favored combina- 
tion, or right-sign (RS) pion, is used for both D° (K s tt~) 
and D° (K s n + ) decays. The squared invariant mass 
of the doubly-Cabibbo-suppressed or wrong-sign (WS) 
two-body combination A"g7r ± (WS) is a linear function of 
M 2 + , DO x and M 2 + _. Three types of final states con- 

tribute in the decay, Cabibbo-favored, doubly-Cabibbo- 
suppressed, and CP eigenstates. The dominant decay 
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FIG. 1: The M(Kg-K + w~) and AM data distributions (points 
with error bars) with fit results overlaid. The dashed lines 
correspond to the background contributions. The vertical 
lines indicate signal ranges used for further analysis. The 
M(Kgiv + iT~) (AM) distribution contains only candidates 
populating the signal AM (M(KgTT + ir~)) range. 



mode is the Cabibbo-favored D° -> K* (892)" tt+, which 
amounts for about 60% of the total branching fraction. 
The second largest contribution is from the intermedi- 
ate CP eigenstate Kgp(770), which is color-suppressed 
compared to K *(892)"7r + . 



VI. DALITZ ANALYSIS 



A simultaneous fit to the resonant substructure of the 
decay to Kg tt + ir~ is performed on the combined D° 
and D° samples to determine the sizes of the various 
contributions. These are compared with previous results 
from other experiments to build confidence in the fit- 
ting technique. Then the fit is applied independently to 
D*-tagged D° and D° samples to measure CP-violation 
asymmetries in the decay amplitudes for each subprocess. 



FIG. 2: Dalitz plot of the decay D°/D° -> K% tt + tt~, where 
the squared invariant masses of the two-body combinations 
J tY§7r ± (RS) and ■K + n are chosen as kinematic quantities. The 
solid line indicates the kinematic boundaries. 



A. Fit to the combined sample 

A binned maximum-likelihood fit to the two- 
dimensional Dalitz-plot distribution with bin widths of 
0.025 GeV 2 /c 4 in both dimensions is performed to deter- 
mine the contributions of the various intermediate reso- 
nances. The likelihood function has the general form 



m = n 



Mi 



(1) 



where 9 are the estimated parameters, I is the number 
of bins, rii is the number of entries in bin i, and is 
the expected number of entries in bin i. The latter are 
obtained using the function 



N[Te(M 2 K o wHBS) ,M* + „-) 
\M(M 2 K „^ (RS) ,M^)\ 2 
(l-T)e(M^ ±(RS) ,M^_) 
\M(M 2 K>Hws) ,M 2 +7T -)\ 2 } 



B { M K%^±(^S)^ M ; 



(2) 



where M. is the complex matrix element of the decay, 
(1 — T) is the fraction of D° candidates with wrongly- 
determined production flavor, called the mistag frac- 
tion, N is the normalization of the number of sig- 
nal events, e {M 2 ( o v ±( RSl yM 2 + ^-) is the relative trigger 
and reconstruction efficiency over the Dalitz plot, and 
B(M^ 0n± , RS y M 2 +n _) is the background distribution. 

The function is evaluated at the bin center to calculate 
the expectation for 



The isobar model [24| is used to describe the matrix 
element M. . The various resonances are modeled by com- 
plex numbers a,je lSj , where j refers to the jth isobar com- 
posed of the amplitude atj and the phase Sj, multiplied 
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with the individual complex matrix element Aj(M,T), 
which depends on the mass M and decay width T of the 
resonance. The phase convention is the same as in the 
CLEO analysis Q and described in Ref. The com- 
plex numbers are added as 

M = a e iSo + aje^AjiM, T) , (3) 

3 

where aoe lS ° represents a possible non-resonant contribu- 
tion. Since aj and Sj are relative amplitudes and phases, 
one resonance can be chosen as the reference. The am- 
plitude and phase of the p(770), being the largest color- 
suppressed mode, are fixed to the values a p (77o) = 1 and 
<5p(770) = 0, respectively. The individual matrix elements, 
Aj(M, T), correspond to normalized Breit-Wigner shapes 
with Blatt-Weisskopf form factor s |25| . A more detailed 
description can be found in Ref. [261 ] . For the intermedi- 
ate resonances p(770) and p(1450) decaying to tt + tt~ , the 
Breit-Wigner shape is replaced by the Gounaris-Sakurai 
description (27| . 

To account for the limited accuracy of the knowledge 
on the masses and widths of the intermediate resonances, 
these parameters can vary within their experimental un- 
certainties, taken from Re f. [ 2q for the /o(980) and 
/o(1370) mesons, and Ref. [29j for the others. This is 
accomplished by means of Gaussian constraints in the 
likelihood function, except for the resonances If* (892) , 
./o(600), and <72> which are unconstrained. Because 
the K*(892) ± is the most prominent resonance, with 
its 60% branching fraction, floating the K*(892) ± mass 
and width in the Dalitz-plot fit improves the fit quality. 
The reason for the unconstrained /o(600) and oi reso- 
nance parameters is the poorly known nature of these 
states. The scalar resonance 02 is introduced to account 
for a structure near lGeV 2 /c 4 in the M 2 , _ distribu- 
tion. A possible explanation for this structure, proposed 
in Ref. [3(|, is the decay /o(980) ^ 77 77 with rescat- 
tcring of to 7r + 7r~, resulting in a distortion of the 
/o(980) -> 7T+ 7r - amplitude for M 2 +7r _ near the rjrj pro- 
duction threshold. The masses and widths of the res- 
onances K*(892) ± , 7^(1430)^ and iY 2 *(1430) ± are re- 
quired to be identical for Cabibbo-favored and doubly- 
Cabibbo-suppressed (DCS) processes. 

Simulated events are used to estimate the relative re- 
construction efficiency over the D° — > K 3 tt + tt~ Dalitz 
plot. The considered decay chain, starting with D* + , is 
simulated by means of the evtgen package [3l|, where 
the three-body decay structure of the D° is generated 
without any intermediate resonances. The generated 
events are passed through the detector simulation and 
reconstructed as data. The simulated detector and trig- 
ger acceptance influence the Dalitz-plot distribution in 
a complicated way. To estimate the efficiency, a binned 
maximum likelihood fit to the Dalitz-plot distribution 
of simulated decays is performed, where a binning of 
0.05 GeV 2 /c 4 in both dimensions is used. An empiric 
function consisting of the sum of a ninth-order multino- 
mial in (M^ 0n±(RS) ) m (M^.) n , where m + n < 9, and 



a Gaussian function, 

e = E a + E x M 2 K a^ ±{KS) + E y M 2 + ^- + E X 2{M^ ±{KS) ) 2 
+ E xv Ml %vHKS) Ml + ^ + E y ,(M 2 + ^) 2 + ... 

(4) 

is employed. The subscripts x and y are abbreviations 
for ^o^±( RS ) and M 2 +7r _, respectively. The Gaussian 

function G(M 2 +n _) models the efficiency enhancement 
at low M 2 + _ values, which is caused by a trigger con- 
figuration that selects track pairs with small opening an- 
gle. The fit projections together with the corresponding 
mass-squared distributions of the three two-body combi- 
nations are shown in Fig. [3] 

The background in the Dalitz-plot distribution re- 
ceives three dominant contributions, combinatorial back- 
ground of purely random particle combinations, misre- 
constructed D° candidates peaking below the D* + signal, 
and combinations of true D candidates with a random 
pion. The first two are estimated in a sample chosen from 
the £>°-mass upper sideband 1.92 < M(K^tt + it~) < 
1.95GeV/c 2 with the same selection requirements used 
in the signal region. The combinations of true D° candi- 
dates with a random pion are directly determined as the 
mistag fraction by the Dalitz-plot fit. 

To estimate the contributions of the individual reso- 
nances to the total decay rate, the fit fractions, 

lKe^A r \ 2 dM 2 K>±{RS) dM 2 + ^ 
r 1 1 Ej ^A 3 \ 2 dM 2 K ^ ±{RS) dM 2 + ^ ' ( } 

are calculated from the fitted amplitudes and phases. 
The statistical uncertainties on the fit fractions are de- 
termined by propagating the uncertainties on the am- 
plitudes and phases. This is done by generating 1000 
random parameter sets of amplitudes and phases accord- 
ing to the full covariance matrix of the fit and taking the 
standard deviation of the distribution of the 1000 calcu- 
lated fit fractions. 

The results of the combined D° and D° Dalitz-plot fit 
for the relative amplitudes and phases of the included 
intermediate resonances can be found in Table HI to- 
gether with the corresponding fit fractions. Table |TT] 
shows the results for the fitted masses and widths of 
the aT*(892) ± , / (600), and cr 2 contributions. The val- 
ues for the K*(892) ± agree with the world-average val- 
ues [29| within 2MeV/c 2 . The mistag fraction obtained 
from the Dalitz-plot fit is 1 - T = (0.98 ± 0.14)%. A 
reduced x 2 of 7387/5082, calculated from the deviations 
between data and fit in each bin, supports the quality 
of our model. The largest discrepancy comes from the 
high statistics corner of the Dalitz plot populated by the 
Cabibbo-favored decays with K*{892) ± resonance. The 
three mass-squared projections are shown in Fig. 2J The 
results for the fit fractions are consistent with the mea- 
surements from previous experiments 0, B [30l [32| . 
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FIG. 3: Squared invariant mass distributions of the simulated 
events used for the determination of the relative reconstruc- 
tion efficiency over the Dalitz plot, together with the corre- 
sponding fit projections. 



B. Measurement of CP-violation asymmetries 

As described in Sec. HI D* tagging is used to measure 
CP- violation effects in the Dalitz decay. Although equal 
numbers of D and D° mesons are produced in the CDF 
II detector, the efficiency for reconstructing soft pions 
from the D* decays causes an instrumental asymmetry 
between the numbers of observed D° and D° decays. 



TABLE I: Combined D° and D° Dalitz-plot-fit results for 
the relative amplitudes and phases of the included intermedi- 
ate resonances, together with the fit fractions calculated from 
them. Due to interference effects between the various reso- 
nances the fit fractions are not constrained to add up exactly 
to 100%. 



Resonance 



Fit fractions 



if*(892) ± 

#5(1430)* 

Jf|(1430) ± 

_R"*(1410) ± 

p(770) 

w(782) 

/o(980) 

/ 2 (1270) 

/o(1370) 

p(1450) 

/o(600) 

oi 

^(892)* (DCS) 
if *(1430) ± (DCS) 
^1(1430)* (DCS) 
Non-Resonant 



1.911 
2.093 
0.986 
1.092 

0.038 
0.476 
1.713 
0.342 
0.709 
1.134 
0.282 
0.137 
0.439 
0.291 
1.797 



±0.012 
± 0.065 
± 0.034 
± 0.069 
1 

± 0.002 
±0.016 
± 0.048 
± 0.021 
± 0.043 
± 0.041 
± 0.023 
± 0.007 
± 0.035 
± 0.034 
± 0.147 



132.1 ±0.7 
54.2 ± 1.9 
308.6 ±2.1 
155.9 ± 2.8 


107.9 ± 2.3 

182.8 ± 1.3 

329.9 ± 1.6 
109.3 ±3.1 

8.7 ±2.7 

201.0 ± 2.9 
16.2 ±9.0 
317.6 ±2.8 

156.1 ±4.9 
213.5 ±6.1 
94.0 ±5.3 



61.80 ±0.31 
6.25 ±0.25 
1.28 ±0.08 
1.07 ±0.10 
18.85 ±0.18 
0.46 ±0.05 
4.91 ±0.19 
1.95 ±0.10 
0.57 ±0.05 
0.41 ±0.04 
7.02 ±0.30 
0.33 ±0.04 
0.32 ±0.03 
0.28 ±0.04 
0.11 ±0.03 
1.64 ± 0.27 



Sum 



107.25 ±0.65 



TABLE II: Combined D° and D° Dalitz-plot-fit results for 
the masses and widths of the J C S'*(892) ± , /o(600), and 02 con- 
tributions. 



Resonance Mass [MeV/c^] 


Natural width [MeV/c^] 


if* (892)* 893.9 ±0.1 
/o(600) 527.3 ±5.2 
cr 2 1150.5 ± 7.7 


51.9 ±0.2 

308.7 ±8.9 

138.8 ± 7.8 



This instrumental asymmetry is mainly due to the tilt of 
COT cells described in Sec. HH which causes positively- 
and negatively-charged particles to hit the cells at differ- 
ent angles. Since only relative differences between the D° 
and D° Dalitz plots are studied, an absolute efficiency 
difference is expected not to bias the observed physics 
asymmetries. However, an instrumental asymmetry de- 
pending on the transverse momentum of the additional 
pion can lead to efficiency discrepancies that vary over 
the Dalitz plot and has to be taken into account to avoid 
biased results. 

Fig. [S] shows the observed asymmetry, 



A = 



N 



D*+ 



N 



D " 



Nr 



(6) 



between the number of D* + (D a ) and D*~ (D°) candi- 
dates as a function of the transverse momentum of the 
pion from the D** decay. The asymmetry is larger at 
low pt("\d*±)- This means that the efficiency for recon- 
structing a D° or a D° may differ over the Dalitz plot. 
The effect is corrected by reweighting the D° Dalitz plot 
according to the deviations between the pt(""d*±) distri- 
butions for positive and negative pion charges found in 
data. 
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FIG. 4: Projections of the Dalitz-plot fit on the individual 
two-body masses, together with the corresponding distribu- 
tions in data. 



Two different parametrization approaches to measure 
CP-violation asymmetries in a simultaneous Dalitz-plot 
fit to the D° and reweighted D° samples are applied. The 
first one corresponds to an independent parametrization 
of the relative amplitudes and phases in the Dalitz-plot 
fits of the D° and D° samples, respectively. Differences 
in the estimated resonance parameters can then be inter- 
preted as CP-violation effects. The second parametriza- 
tion approach is a simultaneous fit to the D° and D° 
samples, where two additional parameters, representing 



FIG. 5: Asymmetry between the numbers of reconstructed 
D* + and D*~ candidates as a function of the soft pion's pr- 



CP-violation amplitudes and phases, are introduced for 
each resonance. 



Independent D° and D° parametrizations 



The fitting procedure described in Sec. lVI Al is repeated 
with separate parametrizations for the amplitudes and 
phases in the D° and D° samples. By performing a si- 
multaneous D° and D fit, common parameters are used 
for the Gaussian-constrained masses and widths of the 
included resonances, the non-resonant contribution, the 
K* (892) ± , /o(600), and 02 masses and widths, as well as 
the mistag fraction. 

To quantify possible CP-violation effects, the fit- 
fraction asymmetries, 



.4 



FF 



FF 



D° 



FF 



D° 



FF, 



FF 



D° 



(7) 



are calculated for each intermediate resonance, where the 
statistical uncertainties are determined by Gaussian un- 
certainties propagated from the statistical uncertainties 
of the individual fit fractions. 

A measure for the overall integrated CP asymmetry is 
given by 



.4 



CP 



\M\ 2 -\M\'' 
\M\ 2 + \M\< 



dM 



if°7r±(RS) 



dM 2 



(8) 



where M. is the matrix element of Eq. ^ for the D° 
decay and M. the one for the D° decay. The statis- 
tical uncertainty on Aqp is determined with the same 
procedure used for the determination of the fit-fraction 
uncertainties. 
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2. CP-violation amplitudes and phases 

Following Ref. 6] , a simultaneous fit to the D° and D° 
samples is performed, where the matrix elements for D° 
and D° read 

M = a e i5 ° + OjeW'+M (l + Aj , 

M = a e i5 ° + £ ajeW'-M I 1 - f J Aj . 
j 

Compared to Eq. Q the additional parameters bj and 
cj>j, representing CP- violation amplitudes and phases, 
are introduced. Again, common parameters are used 
for the Gaussian-constrained masses and widths of the 
included resonances, the non-resonant contribution, the 
K*(892) ± , / (600), and u% masses and widths, as well as 
the mistag fraction. 

VII. SYSTEMATIC UNCERTAINTIES 



in D° and D° data, the Dalitz-plot fits are repeated 
with two independent background samples distinguished 
by the charge of the slow pion in the D*^ decay. The 
systematic uncertainties are calculated as differences be- 
tween the resulting values and the ones from the default 
fits. 



C. Fit model 

The systematic uncertainties originating from the spe- 
cific model used for the Dalitz-plot fit are estimated by 
repeating the fits when one of the resonances fC*(1410) ± , 
/ (1370), cr 2 , X 2 *(1430) ± (DCS), or the non-resonant con- 
tribution is excluded from the model. These are the least 
significant contributions, or in the case of 02, because of 
the controversy over the existence of such a state. For 
each case, the resulting CP-violation quantities are com- 
pared to the values from the default fits and the largest 
deviations are used as modeling systematic uncertainties. 



The systematic uncertainties are categorized into ex- 
perimental and modeling uncertainties. The considered 
experimental sources are efficiency asymmetries varying 
over the Dalitz plot, asymmetries of the background in 
the D° and D° samples, and the applied efficiency distri- 
bution which is estimated by simulated events and may 
not adequately model the composition of trigger configu- 
rations in data. Modeling uncertainties arise from the 
chosen values for the radius parameters in the Blatt- 
Weisskopf form factors and the limited knowledge on the 
complex dynamics of the three-body decay. In this con- 
text, the stability of the determined CP- violation quan- 
tities under variations of the employed Dalitz model is 
tested. The contributions from the various sources to 
the total systematic uncertainties can be found in Ta- 
bles UTTHVU 

A. Efficiency asymmetry 

The reweighting procedure of the D° Dalitz plot ac- 
cording to the deviations between the pt(td*±) distribu- 
tions for positively and negatively charged pions may not 
fully correct for residual small asymmetries between the 
D° and D° efficiency distributions. To estimate the size 
of a systematic effect originating from such an asymme- 
try, the Dalitz-plot fits arc repeated without reweighting 
the D° Dalitz plot. The scale of systematic uncertainties 
are estimated as the differences between the resulting val- 
ues and the ones from the default fits. 



D. Efficiency model 

As described in Sec. [Til different trigger configurations 
with thresholds of 4, 5.5, and 6.5GeV/c on the scalar 
sum of the transverse momenta of the two trigger tracks 
are applied for the online event selection. The simulated 
events, used to determine the reconstruction efficiency 
across the Dalitz plot, are required to pass the different 
trigger configurations in the same proportions as the ac- 
tual data. Any mismodeling of the efficiency is expected 
to cancel between D° and D° candidates, to first order. 
To estimate residual higher-order effects, an efficiency de- 
termined with simulated events satisfying the 6.5 GeV/c 
trigger threshold is used, and the difference with the de- 
fault efficiency is taken as a systematic uncertainty. 

E. Blatt-Weisskopf form factors 

In the default Dalitz-plot fits, the chosen values for the 
radius parameter R in the Blatt-Weisskopf form factors 
are R = 5 (GeV/c)- 1 for the D° and R = 1.5 (GeV/cT 1 
for all intermediate resonances [26j . To estimate the sys- 
tematic uncertainties originating from deviations from 
these values, the Dalitz-plot fits are repeated with both 
values for the radius parameter halved and doubled. The 
systematic uncertainties are then calculated as the largest 
differences between the resulting values and the ones from 
the default fits. 



F. Further checks 

B. Background asymmetry 

To assess the robustness of the fitting procedure, the 
To investigate a possible systematic effect originating following checks were done, each leading to results that 
from different Dalitz-plot distributions of the background agree with the default fits within the expected variations. 
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TABLE III: Systematic uncertainties on the fit-fraction asymmetry Aff for each included intermediate resonance. The contri- 
butions from the efficiency asymmetry, the background asymmetry, the fit model, the trigger efficiency, and the Blatt-Weisskopf 
form factors are described in Sec. IVIII The small values for the p(770) resonance are due to the fixing of its parameters in 
the fit. The total systematic uncertainties are given by adding up the various contributions in quadrature. The corresponding 
statistical uncertainties are listed for comparison. 



Total uncertainties 



A FF [1(T 2 ] 


Efficiency Background Fit model Trigger 


Form factors Systematic Statistical 


K* (892)* 


0.09 


0.13 


0.27 


0.09 


0.23 


0.40 


0.33 


#£(1430)* 


0.9 


0.5 


3.2 


1.1 


1.3 


3.8 


2.4 


#1(1430)* 


0.8 


0.3 


3.9 


0.1 


0.6 


4.1 


4.0 


#*(1410) ± 


0.4 


1.0 


6.2 


0.6 


0.9 


6.4 


5.7 


p(770) 


0.05 


0.00 


0.05 


0.00 


0.04 


0.08 


0.50 


w(782) 


0.6 


0.3 


1.1 


1.0 


2.0 


2.6 


6.0 


/o(980) 


0.9 


1.1 


0.8 


0.1 


0.4 


1.6 


2.2 


/a (1270) 


1.0 


1.4 


1.0 


1.5 


1.8 


3.0 


3.4 


/o(1370) 


0.8 


0.4 


6.4 


0.4 


4.1 


7.7 


4.6 


,9(1450) 


0.6 


0.9 


4.7 


6.2 


1.9 


8.1 


5.2 


/o(600) 


0.2 


0.2 


2.9 


0.2 


2.1 


3.6 


2.7 


(Ji 


0.9 


1.9 


2.5 


0.5 


1.9 


3.8 


7.6 


A'*(892) ± (DCS) 


0.3 


0.6 


1.3 


0.1 


1.5 


2.1 


5.7 


KS (1430)* (DCS) 








5 


8 


5 


10 


11 


Kl (MSO^DCS) 





1 


14 


1 


26 


29 


14 



TABLE IV: Systematic uncertainties on the CP-violation amplitude b for each included intermediate resonance. Further 
explanations can be found in the caption of Table IIIII 

Total uncertainties 



b Efficiency Background Fit model Trigger Form factors Systematic Statistical 



A"*(892) ± 


0.011 


0.001 


0.000 


0.001 


0.002 


0.011 


0.004 


# *(1430)± 


0.022 


0.009 


0.029 


0.009 


0.015 


0.041 


0.028 


#1(1430)* 


0.003 


0.004 


0.022 


0.003 


0.001 


0.023 


0.024 


#*(1410) ± 


0.006 


0.000 


0.015 


0.011 


0.007 


0.021 


0.037 


p(770) 


0.007 


0.003 


0.003 


0.001 


0.001 


0.008 


0.006 


o;(782) 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.002 


/o(980) 


0.002 


0.003 


0.001 


0.000 


0.001 


0.004 


0.005 


/ 2 (1270) 


0.003 


0.008 


0.008 


0.003 


0.005 


0.013 


0.037 


/o(1370) 


0.002 


0.005 


0.008 


0.008 


0.017 


0.021 


0.008 


p(1450) 


0.016 


0.023 


0.107 


0.071 


0.032 


0.135 


0.022 


/o(600) 


0.011 


0.001 


0.018 


0.002 


0.013 


0.025 


0.017 




0.002 


0.002 


0.002 


0.001 


0.001 


0.004 


0.012 


#*(892) ± (DCS) 


0.001 


0.001 


0.001 


0.000 


0.001 


0.002 


0.005 


#5(1430)* (DCS) 


0.003 


0.005 


0.030 


0.008 


0.015 


0.035 


0.024 


# 2 * ( 1430) * (DCS) 


0.005 


0.001 


0.006 


0.002 


0.015 


0.017 


0.029 



To test the effects of the Dalitz-plot binning the 
fits are repeated when varying the bin widths from 
0.025 GeV 2 /c 4 to 0.03GeV 2 /c 4 and 0.05GeV 2 /c 4 . 

Discrepancies between data and the Dalitz-plot fit are 
visible in Fig.|H in particular at the peak of the K* (892) ± 
signal. To verify the results, the fits are repeated when 
excluding the Dalitz-plot regions with the largest discrep- 
ancies between the fit values and data. 

Instead of using Gaussian constraints, the masses and 
decay widths of the included resonances are fixed to 
the values in Ref. [H| for jn(980) and /o(1370), and to 
the world average values [29| for the others, except for 
#*(892) ± , / (600), and o% which are still unconstrained 
parameters in the fit. 

The Dalitz-plot fits are repeated with an alternative 



background model, where the combinatorial background 
and the combinations of true D° candidates with a ran- 
dom pion are modeled by D* + sidebands and the dis- 
tribution of misreconstructed D° candidates is modeled 
by means of an inclusive charm simulated dataset. The 
ratio of combinatorial background and combinations of 
true D a candidates with a random pion is assumed to 
be independent of the D* + mass difference distribution. 
Since this assumption is not completely valid, especially 
close to the kinematic threshold, the method is only cho- 
sen as check of the default background model described 
in Sec. NTM 
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TABLE V: Systematic uncertainties on the CP-violation phase cf> for each included intermediate resonance. Further explanations 
can be found in the caption of Table IIIII 

Total uncertainties 



4> [°] Efficiency Background Fit model Trigger Form factors Systematic Statistical 



K* (892)* 


0.0 


0.1 


1.0 


0.2 


0.8 


1.3 


1.4 


#5(1430)* 


0.0 


0.0 


2.1 


0.2 


0.7 


2.2 


1.7 


#1(1430)* 


0.1 


0.2 


1.0 


0.3 


0.4 


1.1 


1.8 


iT (1410)* 


0.1 


0.0 


0.8 


1.3 


1.6 


2.2 


1.9 


p(770) 


0.1 


0.3 


1.1 


0.2 


0.8 


1.4 


1.5 


a;(782) 


0.3 


0.6 


0.9 


0.1 


0.8 


1.4 


2.2 


/o(980) 


0.2 


0.0 


0.8 


0.0 


0.7 


1.1 


1.3 


/ 2 (1270) 


0.2 


0.2 


1.3 


0.7 


1.4 


2.1 


1.9 


/o(1370) 


1.0 


1.2 


1.9 


1.0 


0.9 


2.8 


1.7 


p(1450) 


0.0 


0.6 


2.0 


3.1 


1.2 


3.9 


1.7 


/o(600) 


0.1 


0.2 


1.3 


0.2 


0.3 


1.4 


1.5 




0.2 


0.3 


0.7 


0.4 


0.7 


1.1 


2.9 


#* (892)* (DCS) 


0.1 


0.5 


1.1 


0.1 


0.1 


1.2 


2.3 


# * (1430)* (DCS) 


0.0 


0.3 


3.5 


1.0 


1.3 


3.9 


4.0 


Kl (1430)* (DCS) 


0.6 


0.5 


2.3 


0.5 


1.7 


3.0 


5.3 



TABLE VI: Systematic uncertainties on the overall integrated CP asymmetry. Further explanations can be found in the caption 
of Table Hm 



Effect 


Uncertainty on Acp [10 z \ 


Efficiency 


0.36 


Background 


0.09 


Fit model 


0.37 


Trigger 


0.05 


Form factors 


0.10 


Total systematic 


0.54 


Statistical 


0.57 



VIII. RESULTS 

All CP-violation quantities are found to be consistent 
with zero. The results for the CP- violation amplitudes 
and phases, defined in Eq. (j9|) and obtained from the 
simultaneous fit to the D° and D° Dalitz plots, are dis- 
played in Table IvTlI The fit-fraction asymmetries for the 
intermediate resonances, defined in Eq. ([7]), are listed in 
Table I Villi The overall time-integrated CP asymmetry, 
defined in Eq. (jH)), is determined to be 

A CP = (-0.05 ± 0.57(stat) ± 0.54(syst))% . (10) 

This value includes the contribution from time-integrated 
CP violation in the mixing of the involved K° mesons. 
We determine this contribution with the method de- 
scribed in Ref. [33] to be —0.07%, much smaller than 
the Aqp measurement uncertainty. 

A. Indirect CP violation 

Following the procedure described in Ref. 16], it is 
possible to disentangle indirect from direct CP-violation 
effects by means of the D° decay time distribution. The 
direct and indirect CP asymmetries are related to the 



TABLE VII: Results of the simultaneous D°-D° Dalitz-plot 
fit for the CP-violation amplitudes, b, and phases, <f>. The 
first uncertainties are statistical and the second systematic. 



Resonance 










b 










<t>[ 


1 






#'(892)* 




+0 


.004 


± 


0.004 


±0 


.011 


-0. 


.8 


± 1 


.4 


± 1. 


.3 


# *(1430)± 




+0 


.044 


± 


0.028 


±0 


.041 


-1. 


.8 


± 1 


.7 


±2. 


.2 


# 2 * (1430)* 




+0 


.018 


± 


0.024 


±0 


.023 


-1. 


1 


± 1 


.8 


± 1. 


.1 


#*(1410)* 




-0 


.010 


± 


0.037 


±0 


.021 


-1. 


6 


± 1 


.9 


±2. 


.2 


p(770) 




-0 


.003 


± 


0.006 


±0 


.008 


-0. 


.5 


± 1 


.5 


± 1. 


.4 


cj(782) 




-0 


.003 


± 


0.002 


±0 


.000 


-1. 


.8 


±2 


.2 


± 1. 


.4 


/o(980) 




-0 


.001 


± 


0.005 


±0 


.004 


-0. 


1 


± 1 


.3 


± 1. 


1 


M1270) 




-0 


.035 


± 


0.037 


±0 


.013 


-2. 





± 1 


.9 


±2. 


1 


/o(1370) 




-0 


.002 


± 


0.008 


±0 


.021 


-0. 


1 


± 1 


.7 


±2. 


.8 


p(1450) 




-0 


.016 


± 


0.022 


±0 


.135 


-1. 


7 


± 1 


.7 


±3. 


.9 


/o(600) 




-0 


.012 


± 


0.017 


±0 


.025 


-0. 


.3 


± 1 


.5 


± 1. 


.4 






-0 


.011 


± 


0.012 


±0 


.004 


-0 


2 


±2 


.9 


± 1. 


1 


K* (892)* (DCS) 


+0 


.001 


± 


0.005 


±0 


.002 


-3. 


.8 


±2 


.3 


± 1. 


.2 


# *(1430)± 


(DCS) +0 


.022 


± 


0.024 


±0 


.035 


-3. 


3 


±4 


.0 


±3. 


.9 


Kl (1430)* 


(DCS) 


-0 


.018 


± 


0.029 


±0 


.017 


+4. 


2 


±5 


.3 


±3. 


.0 



time-integrated asymmetry by 

A C p = Af P + —A^p , (11) 
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TABLE VIII: Fit-fraction asymmetries, Aff, for the included 
intermediate resonances. The first uncertainties are statistical 
and the second systematic. 



Resonance 




A FP [%\ 


^T*(892) ± 




+0.36 ±0.33 ±0.40 


if *(1430)± 




±4.0 ±2.4 ±3.8 


Jf 2 *(1430) ± 




±2.9 ±4.0 ±4.1 


JsT^UlOr 




-2.3 ±5.7 ±6.4 


p(770) 




-0.05 ±0.50 ±0.08 


o;(782) 




-12.6 ±6.0 ±2.6 


/o(980) 




-0.4 ±2.2 ±1.6 


/ 2 (1270) 




-4.0 ±3.4 ±3.0 


/o(1370) 




-0.5 ±4.6 ±7.7 


p(1450) 




-4.1 ±5.2 ±8.1 


/o(600) 




-2.7 ±2.7 ±3.6 






-6.8 ±7.6 ±3.8 


if*(892) ± (DCS) 


+1.0 ±5.7 ±2.1 


^(1430)* 


(DCS) 


+12 ± 11 ± 10 




(DCS) 


-10 ± 14 ±29 



where t is the mean D° lifetime. The mean observed 
decay time (t) is determined from the background sub- 

I 



Since these values are positive by construction, only up- 
per limits are given. The calculation is performed with 
the same method used for the determination of the fit 
fractions, where the 90% and 95% quantiles of resulting 
distributions are used as the corresponding C.L. upper 
limits. To account for systematic uncertainties for each 
resonance, the largest values of all fits with the different 
systematic variations are taken. The resulting 90% and 
95% C.L. on the individual CP asymmetries are listed in 
Table [DD 



IX. MODEL-INDEPENDENT APPROACH 

Following Ref. a mo del- independent search for 

CP violation in the Dalitz-plot distribution of the decay 
D — > Kg 7r + 7r~ is performed by comparing the binned 
Dalitz plots for D° and D° meson decays. No assump- 
tions about the resonant substructure of the decay are 
used. The approach serves as a complementary verifica- 
tion of the results from the Dalitz-plot fits described in 
the previous Sections. However, this method only detects 
the presence of a significant CP- violation effect, without 
allowing a determination of the size of the asymmetries. 

The signed significance of the asymmetry between 
the numbers of D° and D° candidates, (Njja — 



tracted D° decay time distribution. We correct for the 
fraction of nonprompt events that is estimated from 
the D* + impact parameter significance distribution, and 
obtain a mean observed decay time of (t) = (2.28 ± 
0.03)t(_D°). To compare with the recent CDF measure- 
ment of CP- violation asymmetries in D° — > ir + ir~ and 
D° ->■ K+K- decays, A^UD -> h+h') = (-0.01 ± 
0.06(stat) ± 0.04(syst))% [la ], we determine the indirect 
CP asymmetry for the case of no direct CP violation to 
be 

= (-0.02 ± 0.25(stat) ± 0.24(syst))% . (12) 



B. Individual CP-violation asymmetries 

The CLEO experiment also quotes CP-violation quan- 
tities called interference fractions, IF, and individual CP 

asymmetries, Acp, in each subresonance [6J. These are 
ip 

defined as AcPj — pjr-, where 



(13) 



TABLE IX: Upper limits for individual CP-violation asym- 
metries. 

Resonance A C p [%} (90% C.L.) A C p [%] (95% C.L.) 



lC(892p 0.014 0.018 

_R" I *(1430) ± 0.80 1.2 

^5(1430)* 0.45 0.62 
A'*(1410) ± 6.6 8.4 

p(770) 0.038 0.051 

w(782) 0.51 0.66 

/o(980) 0.13 0.17 
/ 2 (1270) 1.6 2.1 
/o(1370) 25 37 
p°(1450) 6.5 8.2 

/o(600) 0.17 0.24 
cr 2 3.1 3.9 
/f*(892) ± (DCS) 1.7 2.3 
A' *(1430) ± (DCS) 22 28 

ifl (1430) ± (DCS) 12 16 



Ngo ) / \J N D o + Nqo , is calculated for each bin and stud- 
ied as a function of the squared K^-k^ and 7r + 7r~ masses. 
In this calculation the number of D° events is normalized 
to the one of D°. Possible CP-violation asymmetries 
would appear as clusters of same-sign discrepancies. The 
sum of the squares of the significance asymmetries in each 



LfEfc(W^ sin(0 fc + ^Ak^AjdM^^dM^^ \ 
(I \M\*dM^ ±{RS) dM* +7T _ + / \M 2 dM^ ±(RS) dM^) 
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bin is expected to follow a \ 2 distribution. The p-value 
can be calculated considering the number of degrees of 
freedom equal to the number of Dalitz-plot bins minus 
one (for the normalization). Furthermore, one expects a 
Gaussian distribution with mean equal to and width of 
1 for the histogram of the asymmetry significance distri- 
bution in case of vanishing CP violation. 

The method is verified in simulation and then applied 
to data. As we test relative differences between D° and 
D° at different places in the Dalitz plot, we normalize D° 
and D° to the same area. With this approach, all asym- 
metries that are uniformly distributed over the Dalitz 
plot completely cancel. However, an efficiency asymme- 
try varying over the Dalitz plot may mimic CP violation. 
As described in Sec. lVIBl this problem is also relevant for 
the Dalitz-plot fits, and the reweighting procedure used 
there is applied here as well. 



2 



CD 

t 1-5 



0.5 - 



(a) 




m ^ ( rs> [GeV 2 /c 4 ] 



H =0.003 ±0.014 
ct = 0.987 ±0.009 




Significance 



FIG. 6: Distribution of (a) the asymmetry significance as a 



function of the squared Kgft 
a histogram. 



and 7r 7r masses and (b) as 



The resulting Dalitz-plot distribution of the asymme- 
try significance between the numbers of D° and D° can- 
didates, with bin widths of 0.025 GeV 2 /c 4 in both di- 
mensions, and the corresponding histogram are shown 
in Fig. [6] The parameters obtained from the fit, fi — 



0.003± 0.014 and a = 0.987± 0.009, are consistent with a 
Gaussian distribution centered at zero with unit variance, 
and the p-value calculated from the asymmetry signifi- 
cance distribution is p — 0.96. The model-independent 
approach confirms that no CP violation is observed be- 
tween the D° and D° decay amplitudes into the Kgir + ir~ 
final state. 



X. CONCLUSION 



A Dalitz-amplitude analysis is employed to study the 
resonant substructure of the D° —> K s ir + n~ three-body 
decay. In performing a full Dalitz-plot fit, the relative 
amplitudes, phases, and fit fractions of the various in- 
termediate resonances are determined. The results are 
compatible and comparable in precision to the measure- 
ments from previous experiments @, H, EH H2| . 

In simultaneous fits to the D° and D° Dalitz plots, 
a CP-violation fit fraction, amplitude, and phase are 
determined for each included intermediate resonance. 
None of these is significantly different from zero. This 
also holds for the overall integrated CP asymmetry, 
A CP = (-0.05±0.57(stat)±0.54(syst))%. A complemen- 
tary model-independent search for localized CP- violation 
differences in relative Dalitz-plot densities between the 
binned D° and D° distributions yields a result consis- 
tent with zero, too. In conclusion, the most precise val- 
ues for the overall integrated CP asymmetry as well as 
the CP-violation fit fractions, amplitudes, and phases are 
reported; no indications for any CP-violation effects in 
D° — s- K s 7T + tt~ decays are found, in agreement with the 
standard model. 
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